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£h ! The Hamiltonian of the radial Coulomb Klein-Gordon and second order Dirac 

^ I equations are shown to possess an infinite symmetric tridiagonal matrix struc- 

' ture on the relativistic Coulomb Sturmian basis. This allows us to give an 

analytic representation for the corresponding Coulomb Green's operators in 
terms of continued fractions. The poles of the Green's matrix reproduce the 
exact relativistic hydrogen spectrum. 
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Ph! I- INTRODUCTION 



In quantum mechanics the knowledge of the Green's operator is equivalent to the com- 



plete knowledge of the system. So, having an analytic basis representation for the Green's 
operator can tremendously simplify the actual calculations. If we know the Green's oper- 
ator only of the asymptotic part of the Hamiltonian we can treat the remaining terms as 
^ 1 perturbations and approximate them by finite matrices. 

In a recent publication, Ref. |J, we have proposed a method for calculating matrix repre- 
sentation of Green's operators. If, in some basis representation, the Hamiltonian possesses an 
infinite symmetric tridiagonal (Jacobi) matrix structure the corresponding Green's operator 
can be given in terms of continued fractions. In Ref. |I| , this theorem were exemplified with 
the Green's operator of the non-relativistic Coulomb and harmonic oscillator Hamiltonian, 
and, in Ref. || , an exactly solvable non-relativistic potential problem were considered which 
provides a smooth transition between the Coulomb and the harmonic oscillator problems. 

The aim of this paper is to extend this result for relativistic Coulomb Green's operators, 
i.e., for the Coulomb Green's operator of the Klein-Gordon and of the second order Dirac 
equations. This later is equivalent to the conventional Dirac equation and seems to have 
several advantages. For details see Ref. [[J and references therein. The Coulomb Sturmian 
matrix elements of the second order Dirac equation has already been obtained by Hostler 
H via evaluating complicated contour integrals. Our derivation, however, is much simpler, 
it relies only on the Jacobi-matrix structure of the Hamiltonian, and the result obtained 
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is also better suited for numerical calculations. In Ref. |J the result appears in terms of 
T and hypergeometric functions, while our procedure results in an easily computable and 
analytically continuable continued fraction. 



II. MATRIX ELEMENTS OF RELATIVISTIC COULOMB GREEN'S 

OPERATORS 

The radial Klein-Gordon and second order Dirac equations for Coulomb interaction are 
given by 

h u \e) = o, (i) 

where 

Hu -( E ) 2 „a | 2aZE | £ u ( u + 1 ) _ (2) 
u \hcj he r dr 2 r 2 

Here [i = mc/h, a = e 2 /hc, m is the mass and Z denotes the charge. For the Klein-Gordon 
case u is given by 



« = -I + yi + /(i + l)-(Za)», (3) 
and in the case of the second order Dirac equation for the different spin states we have 



1 1 //. 1 n2 



The relativistic Coulomb Green's operator is defined as the inverse of the Hamiltonian 



1 u ■ 



H U G U — G U H U — l u , (5) 

where l u denote the unit operator of the radial Hilbert space 7i u - 

In complete analogy with the non-relativistic case we can define the relativistic Coulomb 
Sturmian functions as solutions of the Sturm-Liouville problem 



£+* + ^-y*^)W = 0, (6) 



where rj is a real parameter and n = 0, 1, 2, . . . , oo is the radial quantum number. In 
coordinate space representation they take the form, 



(r \nu; rj) 



(n + 2u + 1)! 



{2r 1 r) u+1 e-^L^ +1 {2r 1 r), (7) 



where L is a Laguerre-polinom. The Coulomb Sturmian functions, together with their 
biorthogonal partner (r \nu; rj) — 1/r • (r \nu\ rj), form a basis: i.e., they are orthogonal, 
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{nu; i] \mu; rf) = {nu; rj \mu; rf) = 5 nm , (8) 
and form a complete set in Ti, u 

oo oo 

\nu; rf) {nu; r]\ = \nu; rf) {nu; r)\ = l u . (9) 

n=0 n=0 

A straightforward calculation yields: 



{nu; rjlmu; rf)= — . 

2r] . 



5 nm (2u + 2n + 2) - 6 nm -iJ (n + l)(n + 2u + 2) 



Snm+iJn (2u + n + 1) 



(10) 



Utilizing this relation and considering Eq. (|^) we can easily calculate the Coulomb Sturmian 
matrix elements of H u , 

£Lnm '■= (nu; rj\ H u \mu; rf) = 

r (2azE . . . AE/Tic) 2 - ri 2 + r] 2 \ 
+5 nm {—f— -2(u + n + l)r) + 2(u + n + l) 1 1 ' — U 

f(E/hc) 2 -ri 2 + ri 2 r -\ (11) 

nm ~ 1 1 2 L V(" + + 2u + 2 )J 

/ (E/Tic) 2 - li 2 + rf r— — -\ 

-<Wi I — — L \/n(n + 2U + 1J1, 

which happens to possess a Jacobi-matrix structure. So, the theorem of Ref. [0 is readily 
applicable here. 

Let us consider the oo x oo Green's matrix 

(G u ) nm = (nu; n\G u \mu; rj), (12) 

and let us denote its rank-iV leading principal submatrix by Cr„ . Then, according to Ref. 

en, 

(G^)Yj 1 = H_ij + 5 jN 5 iN KiN+i (13) 
where JF is a continued fraction 

~K i=N I - 1 = , (14) 

V Oj J Oi +N + o 2+ n + ■ ■ • + o n+N H 

whose coefficients are related to the Jacobi matrix 



This continued fraction convergent for bound-state energies, but, by using the method pre- 
sented in Ref. Q, can be continued analytically to the whole complex energy plane. Simple 
matrix inversion gives now the desired Green's matrix. 

In Table I we demonstrate the numerical accuracy of method by evaluating the ground 
and some highly excited sates of relativistic hydrogen-like atoms, which, in fact, correspond 
to the poles of the Dirac Coulomb Green's matrix. In particular, the zeros of the determinant 
of (in) were located. It should be noted that irrespective of the rank N the zeros should 
provide the exact Dirac results. In Table I we have taken 2x2 matrices. Indeed, the results 
of this method, E c f, agree with the exact one in all cases, practically up to the machine 
accuracy, allowing thus to study the fine structure splitting. 

III. SUMMARY 

In this short note we have presented a practical and easy-to-apply procedure for cal- 
culating the Coulomb Sturmian matrix elements of the Coulomb Green's operator of the 
Klein-Gordon and of the second order Dirac equations. The method is relied only on the 
Jacobi-matrix structure of the corresponding Hamiltonians and results in a continued frac- 
tion which can be continued analytically to the whole complex energy plane. 
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TABLES 





, 

energy levels 


171 

A cf 


TP 




hydrogen 
Z = 1 


IS1/2 


-0.5000066521 


-0.5000066521 


-0.5 


2P l/2 


-0.1250020801 


-0.1250020801 


-0.125 


2P3/2 


-0.1250004160 


-0.1250004160 


-0.125 


50P 1/2 


-0.0002000002 


-0.0002000002 


-0.0002 


50P 3/2 


-0.0002000001 


-0.0002000001 


-0.0002 


uranium 


lSl/2 


-4861.1483347 


-4861.1483347 


-4232 


100D 3/2 


-0.4241695002 


-0.4241695002 


-0.4232 


100D 5/2 


-0.4238303306 


-0.4238303306 


-0.4232 



TABLE I. Energy levels of hydrogen-like atoms in atomic units. E c ^ is the relativistic spectrum 
calculated via continued fraction, E q and Eg are textbook values of the relativistic Dirac and 
the non-relativistic Schrodinger spectrum, respectively. 
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